Revealing multiple gaps in the electronic structure of USb 2 using femtosecond optical 

pulses 
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Ultrafast optical spectroscopy is used to study the antiferromagnetic /-electron system USb2. We 
observe the opening of two charge gaps at low temperatures (< 45 K), arising from renormalization 
of the electronic structure. Analysis of our data indicates that one gap is due to hybridization 
between localized /-electron and conduction electron bands, while band renormalization involving 
magnons leads to the emergence of the second gap. These experiments thus enable us to shed light 
on the complex electronic structure emerging at the Fermi surface in /-electron systems. 
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Exotic phenomena, such as unconventional supercon- 
ductivity, the heavy fermion state, or the elusive "hidden 
order" phase of URu2Si2, can emerge from many-body in- 
teractions in /-electron systems [IHjl • These phenomena, 
governed by strong electronic correlation and complex in- 
teractions between the electronic and bosonic degrees of 
freedom, are often associated with the dual nature (lo- 
calized vs. itinerant) of the /-electrons. The itinerant 
response is usually related to the nature of the Fermi 
surface, where very small band renormalization effects, 
captured theoretically through the electron self-energy, 
often escape experimental detection due to lack of res- 
olution. However, these minute changes in the Fermi 
surface cannot be ignored, as they can dramatically mod- 
ify the properties of /-electron materials, leading to the 
emergence of hidden order parameters, increase of the 
effective mass, superconductivity, and the appearance of 
the heavy fermion state [H-0]. A deeper understanding 
of these phenomena thus depends on a detailed knowl- 
edge of the electronic structure near the Fermi surface in 
/-electron systems. 

Antiferromagnetic USb2 [1,0] is an excellent candidate 
for exploring these issues, as it is a moderately correlated 
electron system with a quasi-2D electronic structure that 
exhibits characteristics of both localized and itinerant 5/ 
electrons T3-llj|. Moreover, previous angle-resolved pho- 
toemission spectroscopy (ARPES) studies indicated that 
bosons could play an important role in band renormal- 
ization la-llTf. These studies produced the first mea- 
surements of the self-energy in 5/ electron systems and 
a model for boson-mediated band renormalization. How- 
ever, they could only explore the region near the center 
of the Brillouin zone (k = 0) and could not specify the 
nature of the boson involved. This directly motivated our 
search for new physical properties (e.g., the opening of a 
gap or changes in the quasiparticle (QP) effective mass) 
arising from boson-mediated many-body interactions, in 
order to understand the complex 5/ band structure in 



USb2 over all /c-space. 

Ultrafast optical spectroscopy (UOS) has been quite 
successful in providing such information, offering insight 
into the physics of different strongly correlated materials 
[3], such as superconductors (SCs) flil - Ell ] and heavy 
fermions (HFs) [22|, |23(. In particular, by measuring the 
temperature (T)-dependent QP dynamics and analyzing 



the data with the Rothwarf- Taylor (RT) model 
one can accurately extract small changes in the electronic 
structure (e.g., the opening of a small gap in the DOS), 
even away from the Brillouin zone center. 

Here, we present the first ultrafast time-resolved dif- 
ferential reflectivity AR(t)/R measurements on USb2, 
taken from room temperature down to 6 K. We reveal 
multiple gaps opening in the density of states (DOS), 
associated with emergent QP states at/near the Fermi 
level. We also observe coherent phonon oscillations, 
previously unobserved in metallic uranium systems, the 
analysis of which not only illustrates the temperature 
evolution of the gap structures, but also demonstrates 
that magnons can be involved in band renormalization. 

In our experiment, high quality USb2 samples, with 
Neel temperatures of Tn ~ 203 K, were grown from Sb 
flux and cleaved along the ab plane prior to measurement. 
The transient photoinduced reflectivity change, AR/R, 
was measured using a Ti: sapphire laser oscillator that 
produced ~55 femtosecond (fs) pulses at a center wave- 
length of 830 nm (~1.5 eV) and a repetition rate of 80 
MHz. The pump and probe beams were cross-polarized, 
with fluences of 0.3 /iJ/cm 2 (~5xl0 -5 photoexcited car- 
riers/unit cell) and 0.02 ^J/cm 2 , respectively. The pump 
beam was incident at a ~10 degree angle to the sample 
normal, with the probe beam directed along the normal. 

Figure Q] shows the measured AR/R signals at (a) low 
and (b) high temperatures. Upon photoexcitation, the 
AR/R signal changes nearly instantaneously due to a rise 
in the temperature of the Fermi surface, as excited car- 
riers rapidly equilibrate via electron-electron scattering 
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FIG. 1. Temperature-dependent AR/R data for USb 2 . The 
solid lines show the extracted non-oscillatory background de- 
cay. The insets show the dynamics at short timescales for 
two temperatures (40 K and 200 K). The arrows indicate the 
corresponding decay processes. 

[22| . This initial change reverses sign when the temper- 
ature crosses ^200 K (~ TV). This is consistent with 
the reconstruction of the Fermi surface at Tn, as pre- 
viously seen in dHvA measurements [ll, 14 1. At longer 
delays, after equilibration of the photoexcited carriers, 
the AR/R signals exhibit a damped ultrafast oscillation 
superimposed on a non-oscillating background decay. 

We will first focus on the non-oscillatory relaxation 
for T < Tjv, since it does not significantly change in 
the paramagnetic state (T > Tjv). Below ~ Tn, this 
relaxation can be fitted with three exponential decays 
convoluted with a Gaussian laser pulse G(t) (FHWM~55 
fs): AR/R = [A s e"*/ Ts + A e e^/ T - + A d e^/ T ^ + A ] <g> 
G(t), (A s > 0, A e < 0, A d < 0). Here, A e and r e 
represent a relaxation process that only appears below 
a critical temperature of ~32 K, and r s and r d are 
the time constants of the initial fast decay and very slow 
relaxation, respectively (see Fig. [T]). 

We find that the relaxation characterized by r d has a 
timescale of at least a few hundred picoseconds (ps) at 
all temperatures, which is likely due to thermal diffusion, 
as in similar measurements on other strongly correlated 
Here, we will focus on the other two de- 
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systems 

cay processes that occur on shorter timescales (t <10 ps). 
We clearly observe from Figure [2] that: (a) r s increases 
continuously with T, and shows a sharp upturn at Tjv 
(nearly diverging there); (b) As T decreases through a 
temperature T* (~ 45 K), A s strongly increases, while 
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FIG. 2. T-dependence of amplitudes Aj and excited quasipar- 
ticle densities n 3 T (j — s,e), corresponding to the relaxation 
times t„ ((a) and (b)) and r e ((c) and (d)), respectively. The 
red solid lines are best fits to the data using the RT model. 
The dashed line in (a) is an amplitude fit using a BCS-like 
T-dependence to obtain the value of A 3 (T = 0) for fitting 

t s (t > t*) m. 

t s decreases significantly (inset to Fig. 2(b)); (c) A e and 
r e , which first appear below T\ increase rapidly as T de- 
creases. Similar temperature-dependent behaviour has 
been observed in superconductors 2^, 21 1 and HF com- 
pounds [22], H3, and was attributed to the opening of 
a gap (A) in the DOS. However, the detailed physical 
mechanism behind the gap opening varies for different 
systems and thus demands careful analysis, as described 
below. 

Quasiparticlc relaxation in a system with a narrow gap, 
such as superconductors and HFs, can be approximated 
by the Rothwarf- Taylor model [24|. Here, the recovery 
process is governed by the decay of electrons with en- 
ergies larger than the gap, via the emission of high fre- 
quency bosons (HFBs) that can subsequently rc-cxcitc 
electron-hole pairs (EHPs). The RT model has been ap- 
plied to many materials [21, 22, 25|, using the equations 



n T (T) = A(0)/A(T) - 1, (1) 
l)- 1 + 2?i r ](A + aTA 4 ), (2) 



where tit(T) is the density of thermally excited quasi- 
particles and a, /3 and S are fitting parameters. In 
these equations, we employ the standard form of tit- 
n T oc (TA)Pe" A / T [13, where the choice of p (0 < p < 1) 
depends on the shape of the DOS. 

We can use the RT model to gain more insight into 
the processes characterized by r s and r e . To model both 
of these processes, we use p = 0.5, which represents 
the DOS with a shape similar to that of BCS super- 
conductors, as previously used for heavy fermion sys- 
tems (Fig. 12"]) [2^. We then fit r e (Fig. Hfc,d)) us- 
ing a T- independent constant gap A e , with a value of 
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11.8 mcV. Modeling r s is more complicated, as the sharp 
change across T* necessitates that we fit the data dif- 
ferently above and below T*. At higher temperatures 
(T > T*), we can fit t s using a BCS-like T-dependence 
for the gap: A s ~ 46.1(1 - T/T N ) - 5 (meV). However, at 
lower temperatures (T < T*), t s can be fit with neither 
a BCS-like T-dependent gap nor a constant gap. There- 
fore, we assume a gap with a simple T-dependent form: 
A* = A*(0)(1 - T/T*y>, where 77 is also a fitting pa- 
rameter. With this assumption, we find that t s (T < T*) 
can be reproduced well using A* ~ 12.8(1 - T/T*) - 05 
(meV) . The excellent agreement between the experimen- 
tal results and the curve fits confirms our initial expec- 
tation of a gap opening. 

The quasi-divergence of r s at the Neel temperature 
suggests a quasiparticle gap opening (measured by A s ) 
due to the onset of magnetic order. In USb2, this mag- 
netic order contains contributions from both spin and or- 
bital polarization. Simultaneously, conditions for a pos- 
sible Fermi surface nesting appear because of parallel 
boundaries exis ting between hole and electron pockets 
in this system [ill [13, 14, 28 1. This is similar to the 
previous observation of a spin density wave gap in the 
itinerant antiferromagnctic actinidc UNiGas [271 ]. 

More interestingly, our analysis indicates that two ad- 
ditional gaps open up at lower temperatures (A* and 
A e ), which we will focus on below. Previous ultrafast 
optical studies on HFs (22| have shown that such small 
gaps are typically due to hybridization between the lo- 
calized 5/ electrons and the conduction electrons. How- 
ever, ARPES studies on USb 2 [HOl have revealed an 
extremely narrow band below the Fermi level and a kink- 
like feature in the band dispersion near the zone center 
(Fig. EJc)) that cannot be explained by such hybridiza- 
tion alone. Ref. [l7| shows that electron-boson medi- 
ated processes also contribute to these effects, where 
the bosons participate in interband electron scattering. 
Therefore, we need to consider both of these potential 
contributions in discussing the origin of A* and A e . 

We have performed a theoretical analysis of the elec- 
tronic structure that incorporates both hybridization be- 
tween / electron bands and conduction bands as well 
as boson-mediated processes. To illustrate this, the cal- 
culated results along the T-X direction are shown in 
Figures [3ja) and (b). In our calculations, the non- 
interacting or bare bands are obtained from local density 
approximation (LDA) calculations, and boson-mediated 
band hybridization is introduced via interband scatter- 
ing [13]. Fig. [3] clearly shows that boson-mediated band 
renormalization can lead to multiple gaps and/or kink- 
like structures away from zone-center at/near the Fermi 
level. More specifically, our calculations produce indi- 
rect charge gaps at the Fermi level with a magnitude of 
~10-15 meV for both hole- and electron-like bands. 

These predicted values agree well with current experi- 
mental findings. Previously, it was not possible to com- 




> 


1 c 






So 

c .„«_ 








c 

LU 

DO 010 

XJ -0.15- 

C 

ho 










r x r x 

Crystal momentum 



FIG. 3. Multiple gaps and kinks revealed by ARPES mea- 
surements on USb2 at 12 K. Panels (a) and (b) show the cal- 
culated bare bands and the renormalized bands when includ- 
ing boson-mediated band renormalization 12911 . resp ectively. 
Panel (c) shows the original ARPES data [la . [171 ] . includ- 
ing the single gap-and-kink structure. Panel (d) shows data 
from (c) reduced with the 2D curvature method [3(j. Mul- 
tiple gap structures are marked with red arrows, with their 
energy scales indicated by red horizontal lines. Blue arrows 
mark kinks in the dispersion due to band renormalization. 

pare the calculation to the original ARPES data over the 
full Brillouin zone 15l-ll7l], since the features away from 
zone center were very hard to discern due to the rapidly 
decreasing signal intensity for high k values. Here, we 
use a recently introduced data reduction method involv- 
ing the 2D Fermi surface curvature [3(| to identify those 
features (Fig.[3ld)), allowing us to show that the complex 
multi-gap structure predicted by theory indeed agrees 
well with the ARPES data. Thus, we propose that boson- 
mediated many-body interactions play a prominent role 
in band renormalization at/near the Fermi surface, and 
contribute to the low-T gap openings observed here. 
However, understanding the nature of the QP states as- 
sociated with the band gaps A* and A e , as well as the 
type of boson involved, requires further evaluation. 

We can gain more insight on these issues by carefully 
considering the oscillations in the AR/R signal. It is gen- 
erally accepted that these terahertz (THz) frequency os- 
cillations, due to coherent optical phonons, are initiated 
either via the displacive excitation of coherent phonons 
(DECP) 3l| or a photoexcitation-induced Raman pro- 
cess [32( in a strongly absorbing material. More quanti- 
tative insight is given by subtracting the non-oscillatory 
background and Fourier transforming the result, reveal- 
ing only one frequency component (Figure LUJa)). This 
allows us to fit the oscillatory signal with the expres- 



4 



sion (AR/ R) osc = Ae~ Tt sin(2irvt + <f), where T and v 
are the damping rate and frequency, respectively (Fig. 
Eta)). The T-dependence of v and T is shown in Figure 0] 
(c,d). It can clearly be seen that v and T depend almost 
linearly on temperature above T* (~ 45 K) but exhibit 
more complicated behaviour for T < T*. 

The T-dependence of the phonon frequency and damp- 
ing rate is typically explained by the anharmonic effect 
33l435|. This effect, leading to a renormalization of the 



phonon energy and damping, usually includes contribu- 
tions from thermal expansion of the lattice (Griineisen 
law) and anharmonic phonon-phonon coupling. We can 
thus model the experimentally measured v and T, includ- 
ing these contributions, using 33-35J 



w(T) = cj + Au^(T) +Ai[l + 2n(u )], (3) 
r(r) = A 2 [l + 2n(wo)], (4) 



where u> = 2nv, n{uj) = [e Hu / kBT - l]" 1 , and the shift 
Aui^ from thermal expansion is given by Aa/^T) = 

T 

-7 f(a c +2a a )dT' 

u>o[e — 1]. Thermal expansion factors a.i 

(i = a,c) are obtained from Ref. |36[. Fig. @] (c,d) 
demonstrates that the above model can explain the In- 
dependent behaviour of v and T above T*. However, it 
fails to capture the phonon softening in v (relative to the 
model prediction) and the wiggle structure in T, which 
both appear below T*. 
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FIG. 4. (a) Extracted oscillation (AR/R) osc and its fit (red 
line) at 60 K. (b) Temperature evolution of the gaps A e and 
Ag compared with the phonon energy (dotted line) for T < 
T* . (c) The T-dependence of the oscillation frequency v and 
(d) the damping rate V. The red solid curves are fits to the 
data using the anharmonic effect model. 

It is well known that in magnetic materials, phonon 
renormalization can occur not only through the anhar- 
monic effect, but also through magnon- phonon coupling, 
which is intrinsically caused by the modulation of the ex- 
change integral J(r) via lattice vibrations [37(. In USb2, 



the magnetic structure changes from paramagnetic to an- 
tifcrromagnctic at Tjv, which could change the magnon- 
phonon coupling and consequently the phonon softening. 
However, we did not observe a clear change in v{T) at 
Tjv (inset to Fig. [He)), suggesting that this effect does 
not play a significant role. 

More insight is obtained by again noting that the large 
phonon softening occurs below T*, which is the same 
temperature where the gap A* opens, suggesting that 
both phenomena are closely related. We first note that 
A* should be spin- related, since it is derived from r s , 
which diverges when the QP gap opens at Tjv. This is 
further supported by the fact that a low energy spin ex- 
citation appears below T* in specific heat measurements 
(sec the supplementary information). These observations 
indicate that magnons participate in the boson-mediated 
band renormalization below T*, leading to the opening 
of A*. This band renormalization increases the effec- 
tive mass of the bands at the Fermi level with respect 
to the bare band structure, increasing their DOS. This 
in turn increases the density of QPs near the Fermi sur- 
face, which enhances the screening of atomic forces, lead- 
ing to the phonon softening seen in our experiments. 
The observed phonon softening thus provides further ev- 
idence for boson (ma<7/T.on)-mediated band renormaliza- 
tion, where we note that the details of possible coupling 
of the spin component to the itinerant electrons remain 
unknown and need to be investigated further. 

At T f (~ 30 K), the gap A e opens, which is associated 
with an additional band renormalization at the Fermi 
level. As with A*, the opening of this gap could influ- 
ence phonon softening. However, we did not observe any 
clear change in v (T) around (Fig. gfc)). This implies 
that any enhancement of the DOS associated with A e is 
not large enough to modify the phonon frequency, indi- 
cating that boson-mediated band renormalization is not 
the origin of A e . Instead, this gap is likely due to the hy- 
bridization of /-electron and conduction electron bands, 
as in other HF materials. This is supported by the fact 
that A e is constant below T', consistent with previous 



findings in other HFs 



In contrast, A*, which devi- 



ates from a constant T dependence, has a more complex 
origin (i.e. boson-mediated band renormalization) . More 
importantly, this reveals that the QP states at/near the 
Fermi level associated with A* and A e are quite different. 

The damping rate T of a coherent phonon with energy 
E p h (= hv) will be strongly enhanced upon a gap (A) 
opening if E p h > A, due to increased coupling between 
the phononic and electronic degrees of freedom. In con- 
trast, T will be unaffected for E p h < A [H, Here, 
there are two gaps, A* and A e , which successively open 
as T decreases. From Fig. Ub), we can see that: (1) for 
T < Tt, E ph is always larger than A e ; (2) for T > T\ 
Eph is greater than A*. Thus, as shown in Fig. Efd), 
the damping rate T below T* should be always greater 
than that resulting from only considering the anharmonic 
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effect, since E p h is always greater than one of the two 
gaps. In addition, since A* gradually increases as T de- 
creases, the coherent phonon energy E ph < A* below 
~ T T . Therefore, the damping associated with A* de- 
creases as T decreases from T* to (~ 30 K). Based on 
these considerations, the dependence of T on T is qual- 
itatively expected to show the behaviour in Fig. 2£d). 
Clearly, T-dependent phonon damping can intuitively re- 
flect the temperature evolution of the gap structure(s) in 
strongly correlated systems. 

In conclusion, we used ultrafast optical spectroscopy 
to shed light on the detailed electronic structure of 
USb2. Temperature-dependent QP relaxation dynam- 
ics revealed the opening of three gaps at different tem- 
peratures, Tt, T*, and T N . The magnitudes of these 
gaps agree well with previous ARPES results and boson- 
mediated band renormalization calculations. More in- 
sight into their properties is given through the T depen- 
dence of the frequency and damping of coherent optical 
phonons. Strong phonon energy renormalization below 
T* also indicates that magnons are the bosons involved 
in band renormalization, which greatly increases the QP 
effective mass at the Fermi level. Overall, these findings 
significantly enhance our understanding of the complex 
emergent states in USb2, as well as in other /-electron 
systems. 
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